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Abstract — Adaptive control technique is adopted to 
synchronize two identical non-autonomous systems with 
unknown parameters in finite time. A virtual unknown 
parameter is introduced in order to avoid the unknown 
parameters from appearing in the controllers and 
parameters update laws. The Duffing equation and a 
gyrostat system are chosen as the numerical examples to 
show the validity of the present method. 
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I. Introduction 

Chaos synchronization has been an interesting research 
area since the pioneer work by Pecora and Carroll [1], 
partially due to its potential applications [2, 3]. So far, most 
of the results on chaos synchronization focus on 
asymptotically stable synchronization, that is, 
synchronization may be achieved with infinite settling time. 
But in practice, one may concern firstly how to stabilize the 
systems within a finite-time interval. Different control 
methods, such as output feedback control and finite time 
observer, were constructed to stabilize nonlinear systems in 
finite time [4-10]. It was demonstrated that the finite-time 
control techniques have better robustness and disturbance 
rejection properties [11-13]. Recently finite-time control 
techniques have been applied in the study of chaos 
synchronization. An approach based upon general results 
borrowed from the robust control theory was used to 
investigate finite-time global chaos synchronization for 
piecewise linear maps [14]. Sliding mode control was an 
effective method to realize finite -time synchronization 
between two chaotic systems with uncertainties [15-19]. 

Controllers, involving the term -ksign(x)\x\ a with 

< a < 1 or e p with 0< P = ql />< 1 , were other valid 
methods often adopted for the studies of finite -time 
synchronization of two chaotic systems [20, 21]. A new 
method, stabilizing subsystem step by step in finite time, 
was proposed to study finite-time synchronization of unified 
chaotic system with uncertain parameters [21], and similar 



problem was further studied based on CLF (control 
Lyapunov function) and Lie derivative [22]. 

The uncertainties concerned in Refs.f 15-22] were mostly 
unmatched parameters, parameter disturbances or external 
disturbances. So far as we know, less attention has been 
paid to the issue of finite-time synchronization of chaotic 
systems with unknown parameters. In this paper, adaptive 
control technique is adopted to synchronize two identical 
non-autonomous chaotic systems with unknown parameters 
in finite time. While the traditional adaptive control method 
is used to design proper controller to realize finite-time 
synchronization, we found that the controllers and 
parameters update laws have to contain unknown 
parameters, which cannot be implemented in practice. In 
order to overcome this difficulty, a known parameter is 
introduced as a virtual unknown parameter. With the help of 
the virtual unknown parameter, an adaptive controller and 
the corresponding parameters update laws are therefore 
designed to synchronize two coupled chaotic systems in 
finite time. The Duffing equation and a gyrostat system are 
chosen as the numerical examples to show the validity of 
the present method. 

II. Finite-time chaos synchronization scheme 

Consider a class of chaotic systems with unknown 
parameters of the form 

x = f p (t,x) + F p (t,x)0 p , (1) 

where t e[0,co) , xeR" is the state vector , 9 p = 
(0 l ,0 2 ,---,0 m ) T eR m is an unknown parameter vector, 
f p (t,x) is a n-dimensional nonlinear function vector and 
F p (t,x) is a nxm nonlinear function matrix. Let Q. x a R" 

be a bounded region containing the whole attractor of drive 
system (1) such that no trajectory of system (1) ever leaves 
it. This assumption is simply based on the bounded property 
of chaotic attractor. Also, let M c/f™ be the set of 
parameter under which system (1) is in a chaotic state. 

Introducing a virtual unknown parameter m+] , Eq.(l) can 
be transformed into 



x = f{t,x) + F(t,x)0 



(2) 



\9\ < d n i = l,2,---,m + l . 



where 9 = (9 l ,9 2 ,---,9 m ,9 m+l ) T eR m+[ is a new parameter 

vector, f(t,x) a n-dimensional nonlinear function vector 

and F(t,x) is a nx(m + Y) nonlinear function matrix. 

Many chaotic systems of the form of system (1) can be 
transformed into system (2). Let us take Lorenz system as 
an example, 
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where 1 ,0 2 ,0 3 are unknown parameters. By introducing a 
virtual unknown parameter 4 =—l, the coefficient of x 2 , 
the Lorenz system can be rewritten as 
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= f(t,x) + F(t,x)0 
If system (2) is chosen as the drive system, the response 
system is constructed as follow, 

y = f(t,y) + F(t,y)9 + u(t), (3) 
where y e R" is the state vector, / (t, y) is a n-dimensional 
nonlinear function vector and F(t,y) is a «x(m + l) 

nonlinear function matrix., 9 e 7? m+1 is a adaptive parameter 
vector and u{t) is a controller. Let c if" be a bounded 

region containing the whole attractor of response system (3) 
without control. 

The object of this paper is to design a proper adaptive 
controller u(t) such that the response system (3) will 

synchronize the drive system (2) in finite time. To this end, 
some hypotheses are needed: 

HI. The nonlinear functions fit,-) and F(t,-) are 
continuous on a bounded closed region [0,oo)xQ and Q 
containing both Q t and Q v . 

H2. The nonlinear matrix function f(t,-) satisfies Lipschitz 
condition in the region [0,co)xQ, 

\\f(t,x)-f(t,y)\\<L\\x-y\\, 
where L is an appropriate positive constant. In this paper, 
I • I denotes matrix or vector norm, defined as 

m n 

INI = (ZZ^) 1 ' 2 fOT the mafriX A = ■ 

H3.The unknown parameter 0— (0 { ,9 2 ,---,0 m ,9 m+i ) T is 
norm bounded. Suppose that 



then 



IMI - d > 

where d = (df +d 2 +■ 



d m+l ) ' are positive constants. 



Define the error variable as e = x-y . Subtracting (3) 
from (2) yields the synchronization error dynamics as 

e = f(t, x) - f{t, y) + F(t, x)0 - F(t, y)0 -u(t) (4) 
Now the finite-time synchronization between systems (2) 
and (3) is transformed into the finite-time stability of error 
dynamics (4), which means that there exists a constant 
T>0 such that 

lim llell = lim llx - y\\ = 0, 

and ||e|| = ||x - y|| = 0, for t > T . Some necessary lemmas are 

introduced in the follows. 

Lemma 1 [21] Assume that a continuous, positive- 
definite function V{t) satisfies the following differential 
inequality, 

V(t) < -cV" (t) , for any t > 0, V(t ) > 0, 
where c>0,0<?7<1 are all constants. Then, for any given 
t Q , V(t) satisfies the following inequality, 
V l ~ n (0 < V ] -« (t )- c(l - ?j)(t -t Q ),t <t< tl , 

and V(t) = for t > with ?,=?„+ — . 

C(l-T]) 

Based on Lemma 1 and motivated by the Theorem 4.10 in 
[24], we establish a lemma for finite-time stable of non- 
autonomous systems. 

Consider the non-autonomous system 

z = g(t,z) (5) 

where g :[0,oo)x£) — >i?" is piecewise continuous in t and 
locally Lipchitz in z on [0,co)xZ) , and D<^R" is a 
domain that contains the origin z = . The origin is an 
equilibrium point for Eq.(5) at t = if g(t,0) = , for 
t>0. 

Lemma 2 Let z = be an equilibrium point for Eq.(5) 
and D(^R" be a domain containing the origin z = . Let 
V ^[0,oo)x D — > R" be a continuously differentiable 
function such that 



a, z <V(t,z)<a 2 z 



— + — g(t, z )^ ~ a 3 \\ z \\ > 
ot oz 



(6) 



(7) 



for any t>0 and zeD , where aj,a 2 ,a 3 ,/> and q are 

positive constants and — < 1 . Then, z = is finite-time 

P 

stable. If the assumptions hold globally, then z = is 
globally finite-time stable. 



Proof: The derivative of V along the trajectories of (5) is 
given by 

z) = — + — g(* , z) < -flj H < . 
ot Sz 

By the inequalities (6), trajectories starting in 
{a 2 \\zf < h} c {V(t,z) < h} , for sufficiently small h , 
remain bounded for all t > t . Inequalities (6) and (7) show 
that V satisfies the differential inequality 

V<-a % \z\ = -a 3 (\\z\\ p y <-a 3 (— )" = -a 3 a 2 "V p , 

a 2 

V~* (t) < (f) - aa p (i _ % - t) , t <t<t { , 

P 



and V(t) = , for t > t x with t 1 =t + 



a';V "(to) 
« 3 (1--) 



By the inequalities (6), 

so z(t) = , for t > f, . Thus, the origin is finite-time stable. 

If the assumptions hold globally, h can be chosen 
arbitrarily large and the foregoing inequality holds for all 

z(f )e*". 

Theorem 1 The coupled system (2) and (3) can be 
synchronized in finite time if the hypotheses H1-H3 hold 
and the following conditions (I)-(II) are satisfied, 
(I) The controller u(t) is chosen as 

k e 

u(t) = V + inr + F (t, x)0 - F(t, y)0 

k 



2k 2 (d 2 +d\\0\\) 



(8) 



where >L,k 2 > . 

(II) The update law of the parameter is 



= F(t,xfe + 



: . , k 2 (A-0) 



(9) 



where A = (d l ,d 2 ,---,d m+l ) T and m+] is the (m + Y) -th 
component of vectors 0. 

For a fixed t , the finite time t x is determined by 



fl=f ° + F F2(?o) ■ 



(10) 



Where V is a positive-definite function satisfying Lemma 2. 
Proof: Choose a Lyapunov function of the form 

1 T 1 

I 

2 



F(e,0) = 



where 9 = 9-9 . In fact, F(e,<9) = -(||e|| +\0\ ) satisfies 

Eq.(6) in Lemma 2 in the case of z = (e, 9) , that is, 

^(|H| 2 + \\&f) < V(e, 9) < \\ef + \\9~f . (1 1) 

The time-derivative of V along the error system (4) is 
V = e T e-(0-0) T S 
= e T (f(t,x)-f(t,y) + F(t,x)0- 

F(t,y)9-u(t))-(9-9f9 

From above hypotheses H1-H3, Eqs.(8) and (9), one can 
obtain 

V = e T (f(t, x) - f(t, y) + F(t, x)0 - F(t, y)9 
-k x e-^.-F(t,x)O + F(t,y)0- 



2k 2 (d 2 +d 0) 



, k 2 jA-0) 



7 )-(0-0) T (F(t, X fe 
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<Z||e|| -A, e H 



2k 2 (d 2 + d jep k 2 {0-0) T {A-0) 



By using k x > L , it yields that 



m+l m+\ 



k 2 (0-0f(0-A) 2k 2 (d 2 +d\ffp 



(12) 



Using the inequalities 

(9-9f(9-A)<2(d 2 +d\9^ 

and 

we have 

v • : k ( HI i "" " ' ' 



= -Allien + ii^ii) < 7ikir + ii^ir (i3) 

The equation (13) shows that the inequality (7) in Lemma 
2 is fit. By Lemma 2 the error system (4) is stabilized at 
origin point e = in finite time. Furthermore the settling 
time t, can be defined by 

V 2 (t ) 2 - 

t l= t + K -f = t + -v 2 (t ). 



That is, the response system (3) synchronizes the drive 
system (2) in finite time. 

Remark 2. If we use the same method in [23] instead of 
the virtual unknown parameter, the controller and parameter 
update laws should be 

u(t) = v+ tt + F ( x ^ - F (y)Q 



e 

ikir 



2k 2 {d 2 + d 
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which contain unknown parameters i (1 < i < rri) of the 
drive system. Theorem 1 indicates that introduction of a 
virtual unknown parameter can avoid the unknown 
parameters from appearing in controllers and parameters 
update laws and realize the finite-time synchronization. 

Remark 3. The magnitude of c/||e|| in the controller 

u(t) will turn to infinity as e — »0 . So, we take <?/(||e|| 2 +s) 

instead of <?/||<?|| 2 in practice with e a sufficient small 
positive constant. Such substitution is also used in [25]. 

III. Numerical simulation 

The Duffing equation and a gyrostat system are selected 
as examples to verify the validity of the proposed method, 
respectively. 

Example 1 The classical Duffing equation with some 
unknown parameters is assumed to be 

x + l x-x + x 3 = 2 coscot , (14) 
where , 2 and CO are unknown parameters. Under some 
parameter values 6> = 0.25 , 6 2 = 0.4 and a> = \, the Duffing 
equation exhibits chaos behavior. Let 
x, = x,x 2 = x , x 3 =cot , and rewrite Eq.(14) in vector form 
f - A f 
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By introducing a virtual unknown parameter 3 = — 1 , the 
coefficient of xf , the Duffing equation can be rewritten as 
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We chose Eq.(15) as the drive system and the response 
system is constructed as 
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(16) 



The controllers u^(t),u 2 (t) and u 3 (t) are defined by 
Eq.(8). In simulation, values of uncertain parameters of the 
drive system are chosen as 9 X = 0.25,6*2 =0.4, 3 =-1 and 
co = \. The bounds of unknown parameters are assumed to 
be d x =0.5,e? 2 =l,d 3 =1 and d 4 =2 . The small positive 
constant is set to be s = 0.0001 , and k } = 1 and k 2 = 0.0001 . 
The initial values are (x, (0), x 2 (0), x 3 (0)) =(1,1,1) , 
(ji (0), y 2 (0), y 3 (0)) = (-2, 3, -1) , (6> (0), 2 (0), 3 (0), «(0)) = 
(0.1,0.2,0.5,0.3) , respectively. From Fig.l, one can see the 
error between the drive and response systems 

||e|| = ^X[ 2 +x\ +X3 converges to zero at about t = 2. The 
time histories of controllers Mj(£),w 2 (f) and m 3 (^) are 
shown in Fig.2. 




Fig.l Error between the drive system (15) and response 
system (16) 




Fig.2 Time histories of controllers in Eq.(16), solid curve 
for u x (t) , dashed curve for u 2 (t) and dotted curve for 
u 3 (t) 

Example 2 A gyrostat system with some unknown 
parameters is assumed to be 



x, = -x 2 x 3 -(flj +a 2 cost)x 2 + 0.4x 3 -0.002(x, + x, 3 ) 
i 2 = x,x 3 -0.4x 3 +a 2 cos?)Xj -0.002(x 2 + x 2 ) (17) 
x 3 = -0.2x 1 +0.2x 2 -0.2x 3 -0.001(x 3 + x 3 ) + a 3 sin? 
where a, , a 2 and a 3 are unknown parameters. Under some 
parameter values, the gyrostat system exhibits chaos 
behavior. Readers refer to Ref.[26] for more details. 
Introduce a virtual unknown parameter a 4 =0.4 , the 

coefficient of x 3 in the first equation in Eq.(17), and rewrite 
it in vector form 
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(18) 



We chose Eq.(17) as the drive system and the response 
system is constructed as 
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(19) 



In simulation, values of uncertain parameters of the drive 
system are chosen as a, = 0.5, a 2 = 3.25 and a, =1.625. 
The bounds of unknown parameters are assumed to be 
d t =l,d 2 =4,d 3 -2 and d 4 = 0.4 . The small positive 

constant is set to be s = 0.0001 , and k t -4 and 
k 2 =0.0001 . The initial values are (x, (0), x 2 (0), x 3 (0)) 
= (1,1,1) , ( Jl (0),j 2 (0),^ 3 (0)) = (-1,-1,-1) , (3,(0), a,(0), 
(a 1 (0), a 2 (0), a 3 (0), a A (0)) = (0. 1, 2, 0.5, 0. 1) , respectively. 
From Fig. 3, one can see the error between the drive and 
response systems ||e|| converges to zero at about t = 2.2 . 
The time histories of controllers u^i), u 2 (t) and u 3 (t) are 
shown in Fig. 4. 

IV. Conclusion 

This work presents a general method for synchronizing 
two identical non-autonomous chaotic systems with 
unknown parameters in finite time. The introduction of a 
virtual unknown parameter can avoid the unknown 
parameters from appearing in controllers and parameters 
update laws. Numerical simulations on the basis of the 



Duffing equation and a gyrostat system show that the 
proposed method is effective. 




Fig. 3 Error between the drive system (18) and response 
system (19) 




Fig.4 Time histories of controllers in Eq.(19), solid curve 
for «j(0 , dashed curve for u 2 (t) and dotted curve for 

" 3 (0 
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